Mathematische Nachrichten, 3 November 2008 



M-functions for closed extensions of adjoint pairs of operators with appli- 
cations to elliptic boundary problems 

B. M. Brown 1 *, G. Grubb 2 ■**, and I. G. Wood 3 *** 

' Cardiff School of Computer Science, Cardiff University, Queen's Buildings, 5 The Parade, Roath, Cardiff CF24 3AA, UK 

2 Department of Mathematical Sciences, University of Copenhagen, Universitetsparken 5, DK-2100 K0benhavn, Denmark 

3 Institute of Mathematics and Physics, Aberystwyth University, Ceredigion SY23 3BZ, Wales UK 

Received XXXX, revised XXXX, accepted XXXX 
Published online XXXX 

Key words Closed extension, M -function, abstract boundary spaces, boundary triplets, elliptic PDEs, pseu- 
dodifferential boundary operators, essential spectrum 
MSC (2000) 35J25, 35J30, 35J55, 35P05, 47A10, 47A11 

Dedicated to the memory of Leonid R. Volevic. 

In this paper, we combine results on extensions of operators with recent results on the relation between the 
A/-function and the spectrum, to examine the spectral behaviour of boundary value problems. M-functions are 
defined for general closed extensions, and associated with realisations of elliptic operators. In particular, we 
consider both ODE and PDE examples where it is possible for the operator to possess spectral points that can 
not be detected by the M -function. 
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1 Introduction 

The extension theory for unbounded operators in Hilbert spaces has been studied since at least 1929 when von 
Neumann discovered the so-called Krein extension. There are many applications of a general extension theory to 
problems generated by both ODE and PDE examples. In the case of symmetric ODEs the book of Naimark [36] 
characterises the extensions of the minimal operator by means of a Cayley transform between the deficiency 
spaces and determines all of these extensions by the imposition of explicit boundary conditions. For PDEs, adjoint 
pairs of second order elliptic operators, their extensions and boundary value problems were studied in the paper 
of Vishik [43] while Grubb [15] showed that all closed extensions of the minimal operator can be characterised 
by nonlocal boundary conditions, building on work of Lions and Magenes [29, 30] (cf. also Hormander [23]) for 
elliptic operators. 

The theory of boundary value spaces (also known as boundary triplets) associated with symmetric operators 
has its origins in the work of Kocubet [24] and Gorbachuk and Gorbachuk [14] with developments from many 
authors, (see [6,25, 27,28, 35, 37, 39,41]). In this context, the theory of the Weyl-M-f unction was developed by 
Derkach and Malamud [9, 10], where spectral properties of the operator were investigated via the M-function 
and Krein-type resolvent formulae were established. For adjoint pairs of abstract operators, boundary triplets 
were introduced by Vainerman [42] and Lyantze and Storozh [31]. Many of the results proved for the symmetric 
case have subsequently been extended to this situation: see, for instance, Malamud and Mogilevski [32] for 
adjoint pairs of operators, and Malamud and Mogilevski [33, 34] for adjoint pairs of linear relations. Amrein 
and Pearson [1] generalised several results from the classical Weyl-m-function for the one-dimensional Sturm - 
Liouville problem to the case of Schrodinger operators, calling them M -functions, in particular they were able to 
show nesting results for families of M-functions on spherical exterior domains in R 3 . For a recent contribution 
with applications to PDEs and characterisation of eigenvalues as poles of an operator valued Weyl- M-function, 
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we refer the reader to [7]. Further recent developments in this area can be found in [3, 8, 12, 38]. There has 
also been extensive work on Dirichlet-to-Neumann maps, also sometimes known as Poincare-Steklov operators, 
especially in the inverse problems literature. These operators have physical meaning, associating, for instance, 
a surface current to an applied voltage and are, in some sense, the natural PDE realisation of the abstract M- 
function which appears in the theory of boundary triplets discussed above. 

Systems of PDEs and even ODEs occur naturally in physical applications (reaction-diffusion equations, Max- 
well systems, Dirac systems, Lame systems) and there is much interest in the spectral properties of operators 
generated by these. In Grubb [17] and Geymonat and Grubb [13] such systems are extensively discussed and 
inter alia points of essential spectrum are characterised by failure of ellipticity of the operator or the boundary 
condition. An alternative abstract approach for block operator matrices has also been developed (see for example 
Atkinson et al. [2]). 

In this paper, we shall combine results obtained by Grubb on extensions of operators, see for example [15], 
with recent results obtained by Brown, Marietta, Naboko and Wood [7] on the relation between the M-function 
and the spectrum, to examine the spectral behaviour of boundary value problems. 1\ /-functions are defined for 
general closed extensions, and associated with realisations of elliptic operators. In particular, we shall consider 
both ODE and PDE examples where it is possible for the operator to possess spectral points that can not be 
detected by the M -function (unlike the classical Sturm-Liouville case). 

In PDE cases, the kernel of the maximal realisation has infinite dimension; then unbounded operators between 
boundary spaces must be allowed, and it is important to choose the representations of the boundary mappings 
in an efficient way. We here rely on the calculus of pseudodifferential operators (^do's), as introduced through 
works of Calderon, Zygmund, Mihlin, Kohn, Nirenberg, Hormander, Seeley and others around the 1960's, as 
well as the calculus of pseudodifferential boundary operators (-0dbo's) introduced by Boutet de Monvel [4,5] and 
applied and extended by Grubb [17]- [21] and others. 

Plan of the paper. Section 2 contains a discussion of the abstract theory. We begin by recalling the universal 
parametrization of closed extensions A established in [15], based on an invertible reference operator Ap, and 
show how it applies to operators A — A (by use of techniques from [17]), giving rise to a Krein resolvent formula 
and characterisations of kernels and ranges, in terms of an abstract boundary operator T x : V\ — * W\. Next, we 
connect this with the boundary triplets theory, as presented in [7] . We first show that for the realisation A b defined 
by a boundary condition Tiu = BTqu (for a special choice of Tq, Ti), the holomorphic operator family Mg(A) 
defined for A 6 £?(Ab) is homeomorphic to the inverse of the holomorphic family T x defined for A e g(Ap), 
when both exist. This takes care of a special class of boundary conditions. The idea is now developed further 
to include general extensions by considering mappings between subspaces, as in [15]. In the present context, 
this replaces the need to work with relations. M-functions are now defined for all closed extensions A, and the 
operator families Mj(A) and T A together describe spectral properties of the operator. 

In Section 3, the ideas are implemented for realisations of elliptic operators on smooth domains 51 in n-space. 
For second-order strongly elliptic operators it is shown in detail how T A and Mr(A), for Neumann-type boundary 
conditions 71 u = C70M, are carried over to mappings L x and Mx(A) between Sobolev spaces over dft. The 
general closed realisations give rise to operator families L\ and Ml x (A) between closed subspaces of L2 
For systems and higher-order operators, the normal elliptic boundary conditions give rise to M- and L-functions 
between products of Sobolev spaces over dfl. 

Section 4 addresses the inverse question: Are the spectral properties fully described by Mj(A) and T A ? The 
answer is in the affirmative for A G q(A) U g(Ap), and this is sufficiently informative in many situations. But 
it is not so in general: We show, both by a PDE and an ODE matrix example, that there exist cases where the 
Af-function is holomorphic across points in the essential spectrum of A (and of Ap). 

The authors thank the referees for careful reading of this paper and useful suggestions for improvements. 

2 Universal parametrization and M -functions 

2.1 A universal parametrization of closed extensions 

It is assumed in this paper that there is given a pair of closed, densely defined operators A m i n and A' min in a 
Hilbert space H (a so-called adjoint pair) such that the adjoint of A' min is an extension of A m - m and the adjoint of 
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A m ; n is an extension of A' min ; we call these adjoints A max resp. A' max . Moreover we assume that there is given 
a closed, densely defined operator Ap lying between A m ; n and A max and having a bounded everywhere defined 
inverse. Thus we get: 

C Ap C A 

max • 

g g(Ap), A max = (A' min )*, 
^min c A* 3 c A' max , G g(Ap), A' max = A^ in ; (2.1) 

here g(B) denotes the resolvent set of B. We call Ap the reference operator. Let M. and M 1 denote the sets of 
operators A lying between A m - m and A max , resp. A' lying between AJ nin and A' max . We write Am for Au when 
A G .M, resp. A'v for A'u when A' G A4'. When [7 is a closed subspace of H, we denote by fjj the orthogonal 
projection of / onto U ; the projection map is denoted pr^. 
Denote also 

ker A max = Z, keiA' max = Z'; (2.2) 

and let 

prg = A^Amax : D(A max ) -» f(^), = I - pr a : £>(A max ) -> Z, 

pr,,, = (A;)- 1 ^ : £>(4nax) - D(AJ), pr c , = J - prp, : D(A' max ) - Z'. (2.3) 

Then pr^ and pr^, resp. prp, and pr^,, are complementary projections defining the direct sum decompositions 

D(A mstx ) = D(A P )+Z, resp. D(A' max ) = D(A* )+Z'. (2.4) 

We also write pr^ u = up, pr^ u = uq, etc. 

The above statements are verified in [15], which also showed the abstract Green's formula 

(Au,v) - (u,A'v) = ((Au) z ,,v c ) - {u c ,(A'v) z ), for u G D(A max ),v G D(A' max ); (2.5) 

and we recall that in that paper, all the closed operators in A4 were characterised by abstract boundary conditions: 

Theorem 2.1 There is a one-to-one correspondence between all closed operators A G M. and all operators 
T : V — > W, where V and W are closed subspaces of Z resp. Z' , and T is closed with domain D(T) dense in 
V. Here T : V — ► W is defined from A by 



D(T) = pr c D(A), V = D(T), W = pr c , D(A*), 

Tu c = {Au) w for u G D(A); (2.6) 

and A is defined from T : V — > W by 

D(A) = {u G D(A max ) | u c G D{T), (Au) w = Tu c }; (2.7) 
it can also be described by 

u G D(A) u = v + z + Ap 1 {Tz + f),veD{A min ),zeD(T)JeZ'eW. (2.8) 

All closed subspaces V C Z and W C Z' and all closed densely defined operators T : V — > W are reached 
in this correspondence. 

When A corresponds to T : V — > W, 

ker A = kcr T, 

ranl = ranT+(J?e W), (2.9) 

orthogonal sum. In particular, A is Fredholm if and only ifT is so, with the same kernel and cokernel. 
If A, hence also T, is injective, the inverse satisfies 

A^ 1 = Ap 1 +T~ 1 pr w , defined on ran A. (2.10) 

The adjoint A* corresponds to T* : W — > V in the analogous way. In particular, in the case where A m i n = 
A' min and Ap is selfadjoint ( called the symmetric case), A is selfadjoint if and only ifV — W and T is selfadjoint. 
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Remark 2.2 The characterisation is related to that of Vishik [43], but differs in an important way: Vishik was 
concerned with normally solvable operators A (those with closed range), and his operators between subspaces of 
Z and Z' map in the opposite direction of those in [15], covering only a subset of them. In contrast, the theory 
in [15] allowed the characterisation of all closed operators in M. 

There are also some results in [15, Section II. 3] on non-closed extensions. 

Now consider the situation where a spectral parameter A G C is subtracted from the operators in Ai. When 
A G g(Ap), we have a similar situation as above: 

A min - A c A - A c A max - A, A' min - A c A* fj - A c A' majc - A, (2.11) 

and we use the notation M. \, -M' x , and 

ker(A max - A) = Z x , kcr(A' max - A) = Z' x ; 
w x = {A fj - X)-\A - A), pr* = (A% - X)-\A' - A), 

pr£ = /-pr£, pr£, = I - pr£, . (2.12) 

Then we have an immediate corollary of Theorem 2.1: 

Corollary 2.3 Let A G g(Ap). There is a 1-1 correspondence between the closed operators A — A in M.\ 
and the closed, densely defined operators T x : V\ — > W\, where V\ and W\ are closed subspaces of Z\ resp. 
Z'y- here 



D(T X )= W $D(A), V X = D(T*), W- X = W *D(A*), 
T x u^ = ((A - \)u) W - x foru G D(A), 

i X G D(T X ), 

In this correspondence, 



D(A) = {u G D(A max ) | «i G D(T X ), ((A - X)u) W - x - T x u x }. (2.13) 



ker(,4 - A) = kcrT A , 

ran(i-A) = t&iiT x + (H Q W- x ) , (2.14) 
orthogonal sum. In particular, if X G q(A), 

(A - A)" 1 = (A l3 - A)" 1 + iy^ffCF*)- 1 ww x ■ (2-15) 
Here \x-*y denotes the injection X ^> Y. 

The formula (2.15) can be regarded as a universal "Krein resolvent formula". It relates the resolvents of 
the arbitrary operator A and the reference operator Ap in a straightforward way, such that information on the 
spectrum of A can be deduced from information on T x . Note also the formulas in (2.14), which give not only 
a correspondence between kernel dimensions and range codimensions, but an identification between kernels and 
cokernels themselves. 

The resolvent of A was studied in [17], from which we extract the following additional information. Define, 
for A G g(Ap), the bounded operators on H: 

E x = A max (A p - A)" 1 = I + X(Ap - A)" 1 , 

^ A = (A m ^-X)A- 1 =1 -XA-\ 

E' x = A' max (Ap - A)- 1 = I + X{A* p ~ A)- 1 = (E X Y, 

F rx = (A' milx - AH^r 1 = I X(Ap)- 1 = (F x )*; (2.16) 
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then E x and F x are inverses of one another, and so are E' x and F' x . In particular, the operators restrict to 
homeomorphisms 

E z : Z — > Z\, F% : Z\ — > Z, 

E§i : Z' ^ Z£, i$ : Z^ ^ Z. (2.17) 

Moreover, for u 6 D(A max ), w £ D(A' max ), 

pr x u = E x pr^ u, pr^ u = prg u — A(A^ — A) -1 pr^ u, 

pr x , v = E' x pr c , v, prg, v = pr^, w — X(A*p - A)" 1 pr ? , u. (2.18) 

This was shown in [17, Sect. 2], in the symmetric case, and the (elementary) proofs extend verbatim to 
the general case. Similar mappings occur frequently in the literature on extensions. The following theorem 
extends [17, Prop. 2.6], to the non-symmetric situation, with practically the same proof: 

Theorem 2.4 For A e Q(Ap), define the operator G x from Z to Z' by 

G x z = -\pr z , E x z, ze Z. (2.19) 

Then 

D{T X ) = E X D(T), V x = E X V, W x = E' X W, 

{T x E x v, E' x w) = (Tv, to) + (G x v, to), for v e D(T), weW. (2.20) 

Proof. The first line in (2.20) follows from (2.13) in view of (2.18). The second line in calculated as follows: 

For u £ D(A), w £W, 

(Tu ( ,w) = (Au,w) = (Au,F' x E ,x w) = (F x Au,E' x w) 

= ((A - \)(Ap)~ l Au, E lX w) = ((A - X)uf3, E' x w) 
= ((A - X)u, E' x w) - ((A - X)u c , E' x w) 

= (T x u x ,E' x w) + (\u c ,E rx w) = (T x E x u c ,E' x w) + (\E x u C: w). 
This shows the equation in (2.20) when we set = v. □ 

Denote by E v the restriction of E x to a mapping from V to V\, with inverse Fy, and let similarly Efy be the 
restriction of E' x to a mapping from W to W\, with inverse F&. Then the second line of (2.20) can be written 

{E$)*T X E V = T + G^ w on D{T) c V, (2.21) 

where 

Gv.yv = P T w G X iv-z ■ (2.22) 
Equivalently, 

T x = (F$)*(T + G VtW )F v on D(T X ) c V x . (2.23) 

Then the Krein resolvent formula (2.15) can be made more explicit as follows: 
Corollary 2.5 When A e q{A) n g(Ap), T x is invertible, and 

(T A)-i = E x {T + Gb tW )- l (E$)*. (2.24) 

Hence 

(A-X)- 1 = {A P -X)- 1 + iv x ^h E X (T + Gy W )~ 1 (E^)* pr w „ . (2.25) 
Proof. (2.24) follows from (2.23) by inversion, and insertion in (2.15) shows (2.25). □ 
Note that G v w depends in a simple way on V and W and is independent of T. 
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2.2 Connections between the universal parametrization and boundary triplets 

The setting for boundary triplets in the non-symmetric case is the following, according to [7] (with reference 
to [31], [34]): A min , A max , A' min and A' max are given as in the beginning of Section 2.1, and there is given a pair 
of Hilbert spaces H, AC and two pairs of "boundary operators" 

rA n /rA k 

: D(A mstK ) -> X, : D(A' max ) -> x, (2.26) 

To J K. \T' J H 

bounded with respect to the graph norm and surjective, such that 

(Au, v) - (u, A'v) = (Titi, T' v) H ~ {T u, r[v) K , all u e D{A max ), v G D(A' max ), (2.27) 

and 

D(A min ) = D(A max ) n kcrTi n kcrTo, D(A' min ) = D{A' max ) n kcrT; n ker^. 
Note that under the assumption of (2. 1), the choice 

H = Z', KL = Z, Ti= pr z , A max , T Q = pr f , T[ = pr z A' max , T' = pr c , (2.28) 

defines in view of (2.4) and (2.5) a boundary triplet. 

Following [7], the boundary triplet is used to define operators As G A4 and A' B , G M 1 for any pair of 
operators B G £(X, H), B' G C(H, K) by 

D(Ab) = ker(ri — BTq), D(A' b ,) =kev{T' 1 - B'T' ). (2.29) 

In order to discuss resolvents, [7] assumes that q{Ab) ^ 0, which means that the situation with existence of 
an invertible Ap as in (2. 1) can be obtained at least after subtraction of a spectral parameter Ao- Thus we can build 
on the full assumption (2.1) from now on. (Theorem 2.1 shows that there is an abundance of different invertible 
operators in M. then.) 

Definition 2.6 For A G g(A B ), the M -function M B (A) is defined by 

Ms (A) : ran(Ti - BT Q ) -> AC, M B (Ti - BT Q )u = T Q u for all u G Z x ; 

and for A G g(A' B ,), the A/-function M' B ,(\) is defined similarly by 

M' B , (A) : ran(Ti - B'T' ) -> TL, M' B , (T' x - B'T' )v = T' v for all v G Z' x . 

It is shown in [7] that when g(A B ) ^ 0, 

{A B )*=A' B ,. (2.30) 

We shall set all this in relation to the universal parametrization, when the boundary triplet is chosen as in 
(2.28). We assume (2.28) from now on. 

Concerning A B , note that B is taken as a bounded operator from Z to Z', and that 

D(A B ) = {ue D(A max ) | {Au) z > = Bu c }, 

by definition. This shows that for the operator T : V — > W that A B corresponds to by Theorem 2.1, 

V = Z, W = Z', T = B. 

Note that (2.30) follows from Theorem 2.1. 

When Z and Z' are finite dimensional, all operators B will be bounded. But in the case where dim Z = 
dimZ' = oo, Theorem 2.1 shows that unbounded T's must be allowed, to cover general extensions. Therefore 
we in the following take B closed, densely defined and possibly unbounded, and define A B by 

D(A B ) = {u G D(A max ) | u c G D(B), {Au) z , = Bu c }. (2.31) 
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Lemma 2.7 ranfT! — BT ) = Z'. In fact, any f G Z' can be written as f = (Ti — BT )v — Yiv for 
v G D(Ap) taken equal to f. 

Proof. Let v run through D(Ap). Then Tgv = pr^ v = G D{B), and Av runs through H = (ran A m i n ) © 
Z', so Tiv — BTqv = (Av)z' runs through Z' . In other words, we can take v = A~Z X f, for any / G Z'. □ 

Lemma 2.8 For any A G q{Ab), Mb(A) is well-defined as a mapping from Z' to Z by 

A/ B (A)(Ti - BT )u = Touforall u G Z x with T u G D(B), 

also when D(B) is a subset of Z. In fact, 

M B (X) = pr c (7 - (A B - A) _1 (^4 max - A)) A" 1 \ Z ,^ H . (2.32) 

Proof. In the defining equation, we can now only allow those u = z x G Z\ for which pr^ z x G D(B). 
We first show: When / = (Ti - Sr )w for v = A^f as in Lemma 2.7, then there is a z G Z\ with 
T z x G such that 

/ = (ri - BT )z x . (2.33) 

For, let x = (As — \)~ 1 (A — X)v, then i; — x G Za. Moreover, since a; G D(Ab), pr^ a; G D(B) in view of 
(2.31), and pr^ v = G D(B), as already noted. Hence 

(r : - sr )(i» - x) = (Ti - sr )« - /; 

so we can take z x = v — x. We conclude: 

{(^ - BT )z x | z x G Z x , pr c z x G £>(B)} = Z' . 

Next, we note that the z x in (2.33) is uniquely determined from /. For, if / = and z x solves (2.33), then 
z x G D(Ab) = D(Ab — A), which is linearly independent from Z\ since A G q(Ab), so z x = 0. 

Thus we can for any / G Z' set Mb(X)J = pr^ z x where z x is the unique solution of (2.33); this defines a 
linear mapping Ms (A) from Z' to Z. The procedures used above to construct Mb (A) are summed up in (2.32). 

□ 

Since B is closed, Mb (A) is closed, hence continuous, as a mapping from Z' to Z. 

For a further analysis of Mb(A), assume A G g{Ap). Then the maps E x , F x etc. in (2.16) are defined. Let 
z x G Z\, and consider the defining equation 

Mb{\){{Az x ) z , -B Wq z x ) = pr c z A ; (2.34) 

where pr^ z x is required to lie in D(B). By (2.17) and (2.18), there is a unique z G Z such that z A = -E^; in 
fact 

z x = E\z = pr x z; z = F^z x = pr c z x , 

so the requirement is that z x G E Z D(B). 

Writing (2.34) in terms of z, and using that Az x = \z x , we find: 

z = pr c z A = AI B (\)((Az x ) z , - Sz) = M B {X){{\z x ) z , - Bz) 
= M B {X){{XE x z) z , - Bz) = M B (X)(-G X - B)z, (2.35) 

cf. (2.19). Here we observe that the operator to the right of Mb (A) equals T x from Section 2. 1 up to homeomor- 
phisms: 

-G x - B = -(G x +T) = -(E'£,)*T X E Z , (2.36) 

by (2.21); here T x is invertible from E^D(B) onto Z' x . We conclude that Mb(A) is the inverse of the operator 
in (2.36). We have shown: 
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Theorem 2.9 When the boundary triplet is chosen as in (2.28) and A £ q(Ab) H g(Ap), -Mg(A) equals the 
inverse of B + G x = T + G x , also equal to the inverse ofT x modulo homeomorphisms: 

-Ms(A)" 1 = B + G x =T + G x = (e£)*T x E z . (2.37) 

In particular, M B (X) has range D(B). 

With this insight we have access to the straightforward resolvent formula (2.25), which implies in this case: 
Corollary 2.10 For A £ g(A B ) H g{Ap), 

(Ab-X)- 1 = {Ap - A)" 1 -\ Zx ^ H E z M B {X){E%)* Wz , % . (2.38) 

Remark 2.11 Other Krein-type resolvent formulae for realisations in the general framework of relations can 
be found in [34, Section 5.2]. 

We also have the direct link between null-spaces and ranges (2.14), when merely A £ g{Ap). 

Corollary 2.12 For any X £ g(Ap), 

kei{A B - A) = E x z kcr(B + G A ), 

ran(A B - A) = (F'Jl)* ran(B + G x ) + ran(A mi „ - A). (2.39) 

For A £ g(Ap), this adds valuable information to the results from [7] on the connection between eigenvalues 
of Ab and poles of Mb (A). 

The analysis moreover implies that Mb(A) and M' B , (A) are adjoints, at least when A £ g(Ap). 

Observe that T + G x (and T x ) is well-defined for all A £ g{Ap), whereas M B (A) is well-defined for all 
A £ g{A B )\ the latter fact is useful for other purposes. In this way, the two operator families complement each 
other, and, together, contain much spectral information. 

It is noteworthy that the widely studied boundary triplets theory leads to an operator family whose elements 
are inverses of elements of the operator family generated by Theorem 2.1 - compare with Remark 2.2 on the 
connection with Vishik's theory. 

2.3 The 7\/-function for arbitrary closed extensions 

The above considerations do not fully use the potential of Corollary 2.3, Theorem 2.4 and Corollary 2.5, which 
allow much more general boundary operators T : V — > W. But, inspired by the result in Theorem 2.9, we can in 
fact establish useful M -functions in all these other cases, namely homeomorphic to the inverses of the operators 
T x that exist for A £ g(A) n g(Ap), and extended to exist for all A £ g(A). 

Theorem 2.13 Let A be an arbitrary closed densely defined operator between A m { n and A maK , and let T : 
V — > W be the corresponding operator according to Theorem 2.1. For any A £ g( A) there is a bounded operator 
Mt(A) :W^>V, depending holomorphically on A £ g(A), such that when A £ g(Ap), —M-^(X) is the inverse 
ofT + Gy w , and is homeomorphic to T x (as defined in Section 2.1). It satisfies 

M x {X){{Az x ) w - Tpr c z x ) = pr c z x , (2.40) 
for all z x £ Z\ such that pr^ z x £ D(T). Its definition extends to all X £ g(A) by the formula 

M x (\) - pr c (i - (A — A)~ 1 (^4 max - A)) A" 1 i w ^ H (2.41) 
In particular, the resolvent formula 

(A - A)" 1 = (Ap - A)" 1 - i Vx ^h E x M A {X){E lX Y pr^. (2.42) 
holds when X £ g(A) fl g{Ap). For all X £ g{Ap), 

kcr(l-A) = E x kcr(T + G X iW ), 

ran(i - A) = ran(T + Gy W ) +HQW~ X . (2.43) 
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Proof. Following the lines of proofs of Lemma 2.7 and 2.8, we define Mt(A) satisfying (2.40) as follows: 
Let / G W. Let v = A- 1 }; then pr c v = G D(T), and 

(Av) w ~ T pr c v = Av = f. (2.44) 

Next, let x = {A — A) -1 (A — X)v, then z x = v — x lies in Z\ and satisfies pr^ z x G D(T) (since pr^ v = and 
pr c x G D(T), cf. (2.6)). This z A satisfies 

(A Z A V--Tpr (Z A = /, (2.45) 

in view of (2.44) and the fact that x e D(A). 

Next, observe that for any vector z x G Za with pr^ z x G -D(r) such that (2.45) holds, / = implies z A = 0, 

since such a z A lies in the two linearly independent spaces D(A — A) and Z\. So there is indeed a mapping from 
/ to pr^ z x solving (2.45), for any / G W, defining Afj(A). It is described by (2.41). The holomorphicity in 

A £ q{A) is seen from this formula. 

The mapping is connected with T x (cf. Corollary 2.3) as follows: 

When A G g(A) n g(Ap), then z = pr c z x = F v z x , and z x = E v z, so the vectors z x with pr c z x G D{T) 
constitute the space E X D{T). Calculating as in (2.35) we then find that 

z = pr c z x = M x {\){{Az x ) w - Tz) - M x {\)((\z x ) w - Tz) 
= M x (\)((\E x z) w - Bz) - M x {\){-G^ w - T)z, 

so Mj(A) is the inverse of — (T + Gy W ) : D(T) — » The remaining statements follow from Corollary 2.3 
and Corollary 2.5. □ 

Note that when Mt(A) is considered in a neighbourhood of a spectral point of ^4 in g(Ap), then we have not 
only information on the possibility of a pole of Mr(A), but an inverse T x , from which ker(^4 — A) and ran(A — A) 
can be read off. 



3 Applications to elliptic partial differential operators 
3.1 Preliminaries 

For elliptic operators A defined over an open subset 57 of W 1 , n > 1, the null-space of the maximal realisation 
is infinite dimensional, so that there is much more freedom of choice of boundary spaces and mappings than in 
ODE cases. It is necessary to allow unbounded operators between boundary spaces to obtain a theory covering the 
well-known cases. Moreover, there is the problem of regularity of domains: For a given realisation A representing 
a boundary condition, it is not always certain that A* represents an analogous boundary condition, but this can 
often be assured if D(A) is known to be contained in the most regular Sobolev space H m (£l), where m is the 
order of A\ this holds when the boundary condition is elliptic. 

The theory of pseudodifferential boundary operators (Boutet de Monvel [4], [5], and e.g. Grubb [20]- [22]) is 
known as an efficient tool in the treatment of boundary value problems on smooth sets (we call it the V'dbo cal- 
culus for short, similarly to the customary use of ^do for pseudodifferential operator). A guiding principle in the 
construction of general theories would therefore be to make it possible to use the i/'dbo calculus in applications to 
concrete operators. The ^'dbo calculus is a theory for genuine operators and their approximate solution operators, 
with many structural refinements; it has not been customary to study relations in this context. We therefore find 
it adequate to interpret the realisations of elliptic operators in terms of the theory based on [15], that characterises 
the elements in Ai by operators, rather than relations. 

Let f2 be a smooth subset of M. n , with C°° boundary <9£1 = S, let m be a positive integer, and let A = 
S| a |<m o-a(x)D a be an m-th order differential operator on f2 with coefficients in C°°(il) and uniformly elliptic 
(i.e., the principal symbol a°(x, £) = J2\ a \=m a a( x )£, a is invertible for all x G f2, all £ G R" \ {0}). The 
maximal and minimal realisations in H = La(^) act like A in the distribution sense, with domains defined by 

D(A max ) = {u G Z 2 (fi) | Au G i 2 (0)}, D(A miB ) = if m (Q), (3.1) 
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and it is well-known (from ellipticity arguments) that 

A* = A' A* = A' n ?1 

min max? max mirp 

where A' max and A' min are the analogous operators for the formal adjoint A' of A. 

The operators belonging to M. resp. M! defined as in Section 2.1 are called the realisations of A resp. A'. 

We denote by H s {CI) the Sobolev space over CI of order s, namely the space of restrictions to CI of the elements 
of H s {W l ), which consists of the distributions u e S' such that (1 + |£| 2 ) s/2 u € L 2 (K n ). (5' = S'(W l ) is 
Schwartz' space of temperate distributions.) By Hq(CI) we denote the subspace of i? s (R") of elements supported 
infi. Fors > — i, s— ^ not integer, this space can be identified with the closure of Cq°(H) in H s (Cl), also denoted 
Hq(CI). Sobolev spaces over E, iP(E), are defined by use of local coordinates. We denote 7jU = (9^m)|e, 
where 9„ is the derivative along the interior normal n at E. Here jj maps 
For j < m there is an extension jj : D(A max ) — > H~i~% (E), cf. e.g. Lions and Magenes [30]. 

Let us briefly recall the relevant elements of the ipdbo calculus. In its general form it treats operators 

(P+ + G K\ C x {Cl) N C°°m N ' 
A=[ : x -► x (3.3) 

\ T Sj C°°(E) M C°°(E) A/ ' 

Here T is a generalized frace operator, going from 51 to E; K is a so-called Poisson operator (called a potential 
operator or coboundary operator in some other texts), going from E to CI; S is a pseudodifferential operator on 
E; and G is an operator on O called a singular Green operator, a non-pseudodifferential term that has to be 
included in order to have adequate composition rules. P is a ?/?do defined on an open set O D ft, and P+ is 
its truncation to CI, defined by P + u = r + Pe + u, where r + restricts D'(Q) to V'(Cl), and e + extends locally 
integrable functions on CI by zero on CI \ CI. P is assumed to satisfy the so-called transmission condition at E; 
this holds for the operators derived from elliptic differential operators that we consider here. There are suitable 
Sobolev space mapping properties in terms of the orders of the entering operators. 
A solvable elliptic PDE problem 

Au = f on CI, Tu = ip on E, (3.4) 

enters in this framework by an operator (where we suppress the index + on A since it acts locally) 

'A\ C°°(Cl) N 

: C°°(Cl) N -» x (3.5) 
T J C°°(E) M ' 

(note that M = and M' > 0), with the inverse 

C°°(Cl) N 

(R K) : x ^C 00 ^)^; R = Q+ + G, (3.6) 

C°°(E) M ' 

where Q + + G solves the problem (3.4) with <p = and K solves the problem (3.4) with / = 0. Here Q is a 
parametrixof A on CI (for example in case A = —A, Q is the convolution with c n \x\ 2 ~ n on W 1 when n > 3), and 
G is the correction term needed to make Q + + G map into the functions satisfying the homogeneous boundary 
condition. 

Besides providing a convenient terminology, the tjjdbo calculus has the advantage that it gives complete com- 
position rules: When A and A' are two systems as in (3.3), the composed operator AA' again has this structure. 
In particular, a composition TK gives a tj}6o on E, and a composition KT gives a singular Green operator. 
Compositions TP + , TG and ST give trace operators, compositions P+K, GK and KS give Poisson operators. 
These are the facts that we shall mainly use in the present paper. Details on the V»dbo calculus are found e.g. 
in [21], [22]. 
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3.2 A typical second-order case 

To give an impression of the theory, we begin by studying in some detail the case of a second-order strongly 
elliptic operator A. This part is divided into five subsections. In the first one we introduce boundary triplets 
for the operator A. In the next three subsections we concentrate on the case of "pure conditions", i.e. when 
T : Z — > Z' . For this case, we show in 3.2.2 how T can be identified with an operator L representing a 
Neumann-type boundary condition. In 3.2.3, we study the corresponding A/-function, proving, among other 
things, a Krein-type resolvent formula. Subsection 3.2.4 takes a closer look at problems with elliptic boundary 
conditions. Finally, in 3.2.5, we consider the general case when T : V — > W and V, W are subspaces of Z and 
Z' , respectively. 

3.2.1 Boundary triplets 

We begin by introducing boundary triplets for the case of a second-order strongly elliptic operator A i.e., with 
Rea°(a;,£) > co|£| 2 for x € fi and £ £ R n (cq > 0), taking SI bounded. Let sq(x) be the (nonvanishing) 
coefficient of — d 2 when A is written in normal and tangential coordinates at a boundary point x, then A has the 
Green's formula 

(Au, v)l 2 (h) - (u, A'v) L2 ^ = (so7iu, 7o«)l 2 (s) - (7o«, s jiv + -4 7ov)l 2 (£) > (3-7) 

for u, v £ -ff 2 (Sl), with a suitable first-order differential operator A' Q over S. We denote s 7i = V\, s 7i = v' x . 
A simple example was explained in [7, Section 7], namely 

A = — A + p(x) ■ grad, with formal adjoint A'v = —Av — div(pv); (3.8) 

where p is an n-vector of functions in C°°(fl). 

We can assume, after addition of a constant to A if necessary, that the Dirichlet problem for A is uniquely 
solvable. 

The Dirichlet realisation A 1 is the operator lying in M. with domain 

D{A 1 ) = D(A max ) n H^n) = H 2 (n) n H^n) 

(the last equality follows by elliptic regularity theory); it has £ g{A 1 ). Let 

Z \( A ) = i u e H s (n) | (A - X)u = 0}, (3.9) 

for s £ K. It is known from [30] that the trace operators 70 and 71, hence also V\, extend by continuity to 
continuous maps 

7o : Z{{A) -^H s -i(X), 7l and v x : Z{{A) H s -i (S), (3.10) 

for all s £ K. When A £ g(A^), let : tp ^ u denote the Poisson operator solving the semi-homogeneous 
Dirichlet problem 

(A - X)u = in SI, 7ou = </>. (3.11) 
It maps continuously 

if s -5(£) -> i/ s (Sl), for all s £ M. (3.12) 

Moreover, it maps i/ s_ 5(S) homeomorphically onto Z^(A) for all s £ M, with 70 acting as an inverse there. 
(We introduce below a special notation for the restricted operator when s = 0, see (3.21).) Analogously, there is 
a Poisson operator Kif solving (3.11) with A — A replaced by A' — A, mapping H s ~? (S) homeomorphically 
onto Z^(A'), with 70 acting as an inverse there. 

Now define the Dirichlet-to-Neumann operators for each A £ g(A 1 ), 
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they are a first-order elliptic pseudodifferential operators over S, continuous and Fredholm from H s 2 (£) to 
i? s -!(£) for all s G R (details e.g. in [16]). 

We shall use the notation for general trace maps j3 and 77: 

P^-.pu^riu, ueZi{A), (3.14) 

when this operator is well-defined. 

Introduce the trace operators V and T' (from [15], where they were called M and M') by 

Tu = Vl u-P^ vilo u, T'u=v' 1 u-P®y^u. (3.15) 

Here V maps D(A max ) continuously onto -ff^(E) and can also be written T = v\A~ x A max , and T' has the 
analogous properties. Moreover, a generalised Green's formula is valid for all u G D(A max ), v G D(A' max ): 

(Au,v) L2{n) - (u,A'v) L2{n) = (ru,7o«)i _i - (7 u,r't;)_i i i, (3.16) 
where (•, -) s .- s denotes the duality pairing between iP(£) and Furthermore, 

{Au,w) L2{n) = (IX 70^)1^1 for all w G Z°(A'). (3.17) 
(Cf. [15, Th. Ill 1.2].) 

To achieve L 2 (£)-dualities in the right-hand side of (3.16), one can choose the norms in iP^E) to be 
induced by suitable isometries from the norm in L2(£). There exists a family of pseudodifferential elliptic 
invertible operators A s of order s £ K on S, symmetric with respect to the duality in £2(2) and with A_ s = A" 1 , 
such that when each H S (E) is provided with the norm for which A s is an isometry from fP(E) onto ^2(2), A s 
also maps isometrically onto £f t_s (E), all t, and 

(A_ s¥ 7,A s ^) s _ s = (^,V)i 2 (S) ! </3,^e£ 2 (S). (3.18) 
Then when we introduce composed operators 

r 1 = Air, ri = Air', r = a_i 7o = v' ; (3.19) 

(3.16) takes the form (2.27) with H = K = L 2 (T,): 

Proposition 3.1 For the adjoint pair A m - ln and ^4^^, (3.19) provides a boundary triplet with TL = /C = 

(Au,v) L2 (q) - (u,A'v) L2 ^ = (Tin, rou) i2 ( S) - (T u,T[v) L2 ^, (3.20) 

holds when u £ D(A max ), v G D(A' max ): 

Such reductions to L2-dualities are made in [7] and [37]. [15] did not make the modification by composition 
with A ± 1 , but worked directly with (3.16). (This was in order to avoid introducing too many operators. Another 
reason was that the Sobolev spaces iP(£) do not have a "preferred norm" when s ^ 0; only the duality (•, -) s ,- s 
should be consistent with the self-duality of £2(2). Moreover, when the realisation A represents an elliptic 
boundary condition, D(A) C H 2 (Vl) and the boundary values are in i 2 (S). — Various homeomorphisms were 
used in [17] for the sake of numerical comparison.) 

In the rest of this section, we use the abbreviation H s for if s (S). We shall keep the formulation with dualities 
in the study of pure Neumann-type boundary conditions, but return to (3.19) in connection with more general 
boundary conditions. 

3.2.2 Interpretation of the boundary conditions 

Consider the set-up of Section 2. 1 with Ap = A 1 , the projection pr^ being denoted pr 7 . The realisation A 1 itself 
of course corresponds to the case V" = VK = {0}in Theorem 2.1. 
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Let A be a closed realisation which corresponds to an operator T with V = Z, W — Z' by Theorem 2. 1 . Note 



that 



Z = Zg(A), Z' = Z$(A'), and for A G g{A y ), Z x = Z&A), Z' x = Z^A 1 ); 



closed subspaces of L,2(£i). 

Denote the restrictions of 7 to mappings from Z\ resp. Z 1 ^ to by ^y Zx resp. j z >_ ; they are homeomor- 
phisms 



7z ,:Z x ^H-i, lz , -.Z^H-i, 



(3.21) 



and their inverses j z resp. j z , coincide with resp. K'^ but have the restricted range space. Their adjoints 
map 

1Zx :H^^Z x , r z ,:H^Z' x . 

When A = 0, the A-indications are left out. 

We shall first interpret A in terms of a boundary condition using the maps with A = 0; this stems from [15]. 
The above homeomorphisms allow "translating" an operator T : Z — > Z' to an operator L : H~2 — > H 2 , as in 
the diagram 



Z 



D(L)= 1Z D(T), 



(3.22) 



Z' 



where the horizontal maps are homeomorphisms. In other words, 

L = (H*z>r l Tl z \ with D(L) = lz D(T) 



(3.23) 



(a closed densely defined operator from H 2 to H 2 ). Hereby we have, when ip = ^ z z £ D(L), tjj = ^ Z 'W £ 

(Tz,w) z > = {T^fz V, l z H)z> = (L<p,ip)i_i. 

Note that D(L) = j D(T) = 70 pr ? D(A) = ^ D(A), since 70 vanishes on D(A 1 ). 
Recall the equation defining T from A. 



(Tpr c u, w) = (Au, w) for u £ D(A), w G Z' . 
In view of (3.17), the right-hand side may be written 

(Au,w) = (Tu, 7ou>) 1 1 = (v\u — P 70)1/1 7ou, 7ow) 1 _ 1 for all w G Z'. 
For the left-hand side we have with L defined above, using that 7^ pr^ u = 70 (u — pi' 7 it) = 7o u > 

(Tpr^u,™) = (L"/zW( u ,lz>w)i_i = (L^ u,^ow)i_i. 
Then, when we write jqw — ip, (3.24) takes the form 

(i7 M,V)i _i = (j/iu-P° ^70",^)! _i for all u G D(A), all ^ G 



(3.24) 



(3.25) 



(3.26) 
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Since ip runs through H 2 , this may be written Lj u = v^u — P® ^70"", or, 

viu = (L + P° 0)I/i )7oU, Jqu G D(L). 

So in fact A represents a Neumann-type boundary condition (3.27). 

Conversely, if we want A to represent a given Neumann-type boundary condition 

v\u = C70U, 

where C is a tpdo over E, we see that L has to be taken to act like 

L = C - P° ,„ . 



(3.27) 



(3.28) 



(3.29) 



Now let us turn to the A-dependent case. Here we consider the families A — A and T x and can proceed in a 
very similar way. When working with the concrete boundary Sobolev spaces we find the advantage that Z and 
Z\ are mapped by 70 to the same space H~ 2 . In fact, 



72 A = IzFz, Jz = lz x E Zl and similarly j z >_ = jz'F£, y z > = 1z<^% 
since, e.g., r yoF x u = -fo(u — XA^u) = 7 u, cf. (2.16). 



Let A G q{A 1 ). In the defining equation for T A , 

(T A pr A u, w) = ((A - \)u, w) for u G D(A), w G Z{ 7 



(3.30) 



(3.31) 



we rewrite the two sides as 



(T x u x ,w) = (L A 7 z A ^, 7^)1,-1, 
((A-X)u,w) = (i^iu-P* 7ow,7ow)i 1 



where L x : H 2 — > iJ 2 is defined by 



note that 



i A = (7*, )" 1 t a 7^ a 1 , £(£ A ) = 7z A i?(r A ); 



D{L X ) = 7 o^(T A ) = i E x D{T) = j D(T) = D(L). 



Then since jz x u x = jqu, the operator A — A represents the boundary condition 
L x l0 u = v x u - P^^lou, j u G £(£ A ) = D(L). 
Moreover, in view of Corollary 2.5, L x is related to T + G x as follows: 



T+G 



(3.32) 



(3.33) 



(3.34) 




where the horizontal maps are homeomorphisms. In view of (3.30), they compose to jz resp. (7^, ) 1 , so 

L x = ( lz ,)~ 1 (T + G x ) lz X =L+ (7lO" 1 G X 7z 1 ' (3-35) 
Since D(A — A) = (3.34) and (3.27) define the same boundary condition, hence 



L X = L + P°-P X on D(L) 



(3.36) 
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Remark 3.2 In particular, it can be inferred (e.g. from the case L = 0) that (j z , ) 1 G x ^ z 1 = P° Vl — P x Vl . 
Note how the operator family L x (in this case where V = Z, W = Z') is written as the sum of a A-independent 
operator L (defining the domain of the realisation) and a A-dependent operator P® Vl — P x g Vl , which is universal 
in the sense that it only depends on A, the set f2, and A. It is useful to observe that since G x is continuous from 
Z to Z' for each A, P2 o Ul — P^ , Vl 15 continuous from H~ 2 to H 2 , hence is of order — 1, in contrast to its two 
individual terms that are elliptic of order 1 (having the same principal symbol). 

This analysis results in the theorem: 

Theorem 3.3 For the second-order strongly elliptic operator A introduced above, let Abe a closed realisation 
with pr^ D(A) dense in Z and pr^; D(A*) dense in Z' . Let T : Z — > Z' be the operator it corresponds to by 
Theorem 2.1. 

(i) When Z and Z 1 are mapped to by 70 and Theorem 2.1 is carried over to the setting based on the 
Green's formula (3.16), A corresponds to a closed, densely defined operator L : — > with domain 
D(L) = 7o-D(A) such that A represents the boundary condition 

v x u = C70U, where C = L + P® QMl ■ (3.37) 

Here L is defined from T by (3.23). 

(ii) For any A G q{A 1 ), A — A corresponds similarly to 

L X = L + P° o>Ui - P^ ui : R-\ -> Hi, (3.38) 

with domain D(L X ) = D(L). A — A is the realisation of A — A determined by the boundary condition (3.37), 
where C may also be written C = L x + P X Q Vl . Here, when A — A corresponds to T A by Corollary 2.3, L x is 
defined from T x by (3.32). Moreover, (3.35) holds. 

(iii) Furthermore, for any A 6 q{A~ i ), 

ker(2- A) = K x kevL x , 

T&n(A - A) = i* ZL ran L x + ran(A min - A). (3.39) 
Proof. All has been accounted for above except point (iii), but this follows immediately from (2.14). □ 

3.2.3 The M-function 

We define an ill-function in this representation, by use of Lemma 2.8. We have from (2.32) for A € q(A): 

M A (\) = pr c ((J - (A - A) _1 (^l max - A))^ 1 i z ,^ H :Z'^Z. (3.40) 

We know from Theorem 2.9 that the M-function should coincide with minus the inverse of the operator induced 
from T x in (2.36), when A € q(A) D q{A 1 ). So, applying the trace maps in (3.21) in a similar way as we did for 
T, we get 

Ml (A) = 7zM^(A) 7 *, = lz pr c ((J - (A- A) _1 (A max - A)) A; 1 \ Z ,^ H 7*, 

= 70(7 - {A- X)~ 1 (A max - A))^" 1 i z ^H lz ,. (3.41) 

Here, when A E q(A) n g{A~ f ), — M^(A) is the inverse of the operator translated from T + G x , namely, in view 
of (3.35)-(3.36), 

Ml (A) = —{L + P° o>I/i - P^J- 1 = -(L x )-\ (3.42) 

bounded from Hi to H~*, and then it has range D(L). Moreover, it has the extension to A £ q(A) given in 
(3.41), a holomorphic family of bounded operators from H 2 to H~? . 

Note that when A G q{A) l~l q{A 1 ), L x is surjective onto H 2 . These considerations lead to: 
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Theorem 3.4 For the realisation considered in Theorem 3.3, there is an M -function defined by 

Ml (A) = 70 (I - (A - A) _1 (A max - X^A' 1 i z ^ Hl * z ,, (3.43) 

a holomorphic family of bounded operators from to H~^. For X G g{A) (~l g{A 1 ), it satisfies 

M L {X) = -(L + P 7 ° 0)71 - P^r 1 = -(L x )-\ (3.44) 

There is the following Krein resolvent formula, valid for all X G q{A) D £?(Ay): 

(A - A) - 1 = ( A, - A) - 1 - i ZA ^ l z l Ml ( A) ( 7 £, ) " 1 pr z , (3 .45) 

= (A 1 -X)- 1 -K x M L (\)(K! y x )*. 

Proof. The definition of Ml is accounted for above. The first line in the Krein formula follows from (2.42) 
in the case V = Z,W = Z', by the calculation 

E X Z M X (X)(E> Z \)* = E x ^ML{XW z ,)-\E'lr = T^M^A)^)" 1 , 

using (3.41) and (3.30). The second line follows since iz x ^HJz^ = K x : H~z — > H, (7^, ) _1 pr z , = 
(iz£-H 7zD* = ( K f)* -H^Hi (recall that H = L 2 (fl)). □ 

Note that with the notation (3.14), L x = -P^^-o-ro and M ^( A ) = P i-c 7o , 7 o' cf < 3 ' 29 )- 



3.2.4 Elliptic boundary conditions 

Further information can be obtained in elliptic cases. Using the Sobolev space mapping properties of jz and its 
inverse, one finds since D{A 1 ) c H 2 {VL) that D(A) C H 2 (fl) if and only if D(L) c Hi. If L is given as an 
arbitrary ^do of order 1, it is not in general bounded from iT~3 to H^, but has a subset as domain. 

Ellipticity of the boundary value problem defined from (3.28), (3.29) (the Shapiro-Lopatinskii condition) holds 
precisely when L acts like an elliptic V'do of order 1. Then ranL C H2 implies D(L) C H?, and the graph- 
norm on D(L) is equivalent with the -norm. Let q(A)P\q(A~ ( ) ^ (a reasonable hypothesis in our discussion) 
and let Ao G g(A) (~l q{A~). Then L x ° has kernel and cokernel {0} (since T x ° has so, by Theorem 2.1) and is 
likewise elliptic of order 1, so the inverse — Ml(Xo) is an elliptic ipdo of order — 1; it is defined on all of Hi. It 
maps H 3 onto Hi (since the range of the operator from any H s to H s+1 has codimension 0). It follows that 

D(L) = D(L Xa ) = Hi. (3.46) 

Then moreover, D(A) C H 2 (fl) (and D(A) is the largest subset of D(A max ) where (3.28) holds). 

In this case, the adjoint of (L A °) _1 (as a bounded operator from H 2 to H ~ 2) is ((L A °)*) _1 , so the adjoint of 

11 3 ~ 

L (as an unbounded closed densely defined operator from H~2 to H 3) is L* with domain ffJ. Then A* is the 
realisation of A' defined by the elliptic boundary condition 

i/ lU = (L* + P^jou, (3.47) 

and D(A*) C _ff 2 (il). Note that if A = j4_' and L is elliptic of order 1 and symmetric, A is selfadjoint. 

Before including these observations in a theorem we shall show that also pv Zx , pr z / and their adjoints \z x ^h, 

l Z'- ->fl^ belong to the ipdbo calculus; this will allow a discussion of Ml (A) for A G g(A) \ g{A 1 ). 

Proposition 3.5 Let X G g{A~ f ). Then pr Z '_ acts as the singular Green operator 

pT zk =I-(A max -X)R((A'-X)(A-X),j ,j 1 )(A' [aax -X), (3.48) 
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where R((A' — A) (A — A), 70, 71) : g u is the solution operator for the problem 

{A 1 -X)(A — X)u = g, lo u = jiu = 0. (3.49) 

Similarly, 

V v Zx =/-K iax -A)i?((A-A)(A'-A),7o,7i)(Ana X -A). (3.50) 

Moreover, iz x —>H acts like the adjoint of the operator in (3.50) in the ijjdbo calculus, andizt QCtS like the 
adjoint of the operator in (3.48). 

Proof. When / € H, pr z > f is the second component of / in the decomposition 

H = ran(A min - A) © Z'y 

Write / = fi + fi according to this decomposition, and note that 

(A' - A)/ = (A' - A)/i- (3.51) 

Observe moreover, that v = (A 1 — \)~ 1 fi is the unique element in D(A m [ n ) = Hq(£1) such that (A — X)v = f\. 
In view of (3.51), v moreover solves 

(A'- \){A- X)v = (A'- A)/, 7o i. = 7 lW = 0, (3.52) 

and this solution is unique since (A' — A) (A — A) is formally selfadjoint strongly elliptic with positive minimal 
realisation, hence has a positive Friedrichs extension, representing its Dirichlet problem (3.49). Thus v is uniquely 
determined as 

v = R((A> - X)(A - A), 7o ,7i)(A' - A)/, 

and f-2 is given by the formula (3.48). It should be noted that the operator in (3.48) has a good meaning on L-2,(Q) 
in the ?/>dbo calculus; first A 1 — A maps £2(0) continuously into H~ 2 (£l), then R((A' — A) (A — A), 70, 71) maps 
H~ 2 (£l) homeomorphically onto Hq(CI), as is known for Dirichlet problems for positive operators, and finally 
A — X maps Hq(CI) continuously into L,2(£l). (3.48) defines a singular Green operator since the ipdo part of the 
second term cancels out with /. 

The formula (3.50) follows by interchanging the roles of A — A and A' — A. 

Finally, as a technical point taken care of in [20], the operator in (3.48), although the factor to the right is of 
order 2, is of class since R((A' — X) (A — A), 70, 71) is of class —2. Then it does have an adjoint in the ijjdbo 
calculus, so the assertion follows since the adjoint of pr^ : H — > Z\ is iz\ >h ■ d 

Theorem 3.6 For the operators considered in Theorems 3.3 and 3.4, one has: 

(i) When L acts like an elliptic ipdo of order 1, then D(L) C H 2 and D(A) C H 2 (fl), and when X € 
g(A) fl g{A 1 ), all the operators entering in the formulas belong to the ijjdbo calculus, and Ml[X) is elliptic of 
order -1. Moreover, if q(A) n q{A 7 ) ^ 0, D(L) = ifi. 

(ii) When the statements in (i) hold, we have moreover that A* is the realisation of A' determined by the 
elliptic boundary condition (3.47), and D(L*) = Hz, D(A*) C H 2 (fl). In particular, if A — A' and L is 
symmetric and elliptic of order 1, A is selfadjoint with D(A) C H 2 (Q). 

(iii) When L is elliptic of order 1 with D(L) = Hi, Ml(X) is an elliptic ipdo of order — 1 for all A G q(A). 

Proof. The assertions in (i) and (ii) were shown above. For (iii), we use the formula (3.43) derived in (3.41) 
from (2.32). When L is elliptic with D(L) = H*, the resolvent (A — A) -1 belongs to the ?/>dbo calculus and 
maps Li{H) into H 2 (£l); then the composition rules in the f/;dbo calculus imply that Ml(X) is a ipdo over E. 
Since (J — (A — A)~ 1 (A max — A))^" 1 maps £2^) continuously into H 2 {£1), Ml(X) maps continuously 
into H i, hence is of order — 1, for all A € q(A). Since the principal symbol of Mj(A) (in the ipdho calculus) is 
independent of A, also the principal symbol of the ^do Ml(X) is independent of A; it equals the inverse of the 
principal symbol of L and is therefore elliptic. □ 
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Remark 3.7 The Krein resolvent formula (3.45) can be compared with the standard resolvent formula 
(A — A) -1 = Q\.+ + G\ (cf. Seeley [40] and e.g. Grubb [21]), where Q\ is a parametrix (an approximate 
inverse) of A — A on a neighborhood Q, of 0, Qa,+ is its truncation to ft, and G\ is a singular Green operator 
adapted to the boundary condition (as in (3.6)). This has been used e.g. to show asymptotic kernel and trace 
expansions (also for the associated heat operator). Formula (3.45) gives more direct eigenvalue information. 

Remark 3.8 The operators L and Ml (A) can be pseudodifferential also in non-elliptic cases. A striking case 
is where L = as an operator from H~i to H 2 , so that A represents the boundary condition 

viu = K,^ou. (3.53) 

This is Krein's "soft extension" [26], namely the realisation A m with domain 

D{A M ) = D(A min )+Z, (3.54) 

as described in [15]. Its domain is not contained in any i? s (SI) with s > 0. Here Ml (A) = — (P° Q Vl — P^ Q Ul )~ l 
when this inverse exists; cf. (3.44). Spectral properties of Am are worked out in [19]. 



3.2.5 The general case 

The previous subsections cover the cases where T goes from Z to Z' in Theorem 2.1 (called "pure conditions" 
in [15]). The family M. moreover contains the operators corresponding to closed, densely defined operators 
T : V — > W with arbitrary closed subspaces V C Z, W C Z'. In [15], it was shown that these correspond to 
operators L : X — ► Y*, X = -foV and Y = 70 W (subspaces of H~ 2 ); then A represents the boundary condition 



7oit € X, (L70M, ip)y*,Y = (Tu,ip)i 1, all?/; e Y. 



(3.55) 



The equation is also written Lj u = Tu\y (restriction as a functional on Y), 

Let us show how this looks when we use the modified trace operators in (3.19) mapping the maximal domains 
to L 2 (S). Setting 



X 1 =T V = A_iX,Y 1 =T W = A_xY, L\ = (T, 



o.w. 



^TT-y, withD(L 1 ) = r D(T), 

(3.56) 



where T y resp. T w denote the restrictions of r as mappings from V to X\ resp. from W to Y\, we have the 
following diagram: 



(3.57) 




We 



find as in (3.25)-(3.27) that the statements pr^ u G D(T), Tu^ = (Au)w< carry over to the statements 



T M E D(Li) C Xi, LxTqu = pr Yi Tiu 



(3.58) 

this is then the boundary condition represented by A. Note that the condition Tou G X\ enters as an important 
part of the boundary condition, compensating for the fact that L\ acts between smaller spaces than in the case 
V = Z, W = Z'. 

Since pr Yi Tiu = pr Yl A 1 {ym — P® , Vl lo u )> (3-58) may also be written in terms of the standard traces ■you 
and v\u, as 

A_i 7o u G D(L X ), pr n Aivtu = {L x + pr Yi AiP° o>Vi Ai i Xl ^L 2 )A„iy Q u. (3.59) 
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Such formulations can likewise be pursued for A — A, allowing an extension of Theorems 3.3 and 3.4. Let 

A £ p(Ay). Now r ,y A = T 0)V F^ : V\ ^ X u and r ,^ = T 0tW F$ : ^ Y\- The defining equation for 



T x is: 

{T x pr x u,w) = ((A-X)u,w) foru G D(A),w G W- x . (3.60) 

We define 

L i = (Kw-X'T^ok' withD(L^) = T oyx D(T x ) = T oy D(T) = D(Li), (3.61) 
and then rewrite the two sides as follows, denoting T x = A i (yx — P X Q ^70): 

(t a u^,w) = (L^r 0) v x ^,r , 1 y x M)z, 2 (E) = (L x r u,r w) Yl , 
((A - x)u, w) = (r x u, r w) L2(s) = (pr Yi r x u, r w) Yl . 

When w runs through W^, Toiw runs through Y\, so we see that A — A represents the boundary condition 

T u G L>(Li) C Xi, LJT u = pr Yl T x u. (3.62) 
It can also be written in terms of the standard trace maps as 

A_i 7o u G D(Za) C Xi, L^A_i 7o u = pr yi Ai (vm - P x ^ Q u), (3.63) 
where the equation can be rewritten as 

pr n Ai vm = (L x + pr Yi AiP^Ai i Xl ^z, 2 )A_ 1 7o m. (3.64) 

Since D(2 - A) = L>(2), we have in view of (3.59): 

L X = L 1+ P r Yi Ai(i* iB1 - i^,^)Ai ijf^L,, when A G f>(A 7 ). (3.65) 
All this leads to: 

Theorem 3.9 For f/je second-order strongly elliptic operator A introduced above, let the closed realisation A 
correspond to T : V — » W by Theorem 2.1. 

(i) Define X\, Y\ and L\ by (3.56); f/ien A represents the boundary condition (3.58), more explicitly written 
as (3.59). 

(ii) For any A G q{A 1 ), A — A corresponds to the operator L x : X\ — > Yi defined in (3.61), and represents 
the boundary condition (3.62), also written as in (3.63), (3.64). Here h\ and L\ are related by (3.65), and 

ker(2- A) = K x Ai kcrL x , 

ran(i - A) = T^ w% ranZ^ + (H G W x ). (3.66) 

Proof. In view of the preparations before the theorem, it remains to account for (3.66), which follows by 
application of the various transformation maps to (2.43). □ 

The M -function in this set-up is defined from formula (2.41), in a similar way as in (3.41): 

M Ll (X) = r oy M A (X)T* >w = A_ i7v M A (X)T* tW 

= A_ ijv pr c ((7 -{A- A)- 1 (-4max - A))^ 1 i w ^ H T* ow 

= A_i 7o (J - (A - A) _1 (j4 max - A)) A; 1 i w ^ H T* o w . (3.67) 
Here, when A E q(A) fl q(A 7 ), —Ml 1 (A) is the inverse of L\. We therefore have the following: 
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Theorem 3.10 Once again, let A be the second-order strongly elliptic operator introduced above and let the 
closed realisation A correspond to T : V — > W by Theorem 2.1. 
The M -function in this setting is 

M Ll {\) = A_i 7o (/ - {A-X)-\A max - A)) A; 1 i w ^ H T^ w , (3.68) 

a family of bounded operators from Y\ to X\, depending holomorphically on A € q{A). For A £ q(A) n q(A 1 ),it 
satisfies 

M Ll (X) = -(ii + P r yi Ai (P° iVl - i^.JAi ix^Lj" 1 . (3.69) 
The following resolvent formula holds for all A S q[A) fl g(Ay): 

(I - A)- 1 = (A 7 - A)- 1 - i v ^ H T-^M^WiT*^)- 1 pr Wj . (3.70) 
= (A - A)" 1 - A- 7 A Xi M Ll (A)(< yi )*, 
where Xi : X\ — > 7? acta Z/£e the composition of Ki : X\ — > A, i i and : — > i7. 

Proof. It remains to show (3.70), which follows from (2.42). □ 

The analysis of the realisations corresponding to operators T : Z — > Z' covers all the most frequently studied 
boundary conditions for second-order scalar elliptic operators, whereas the cases where T acts between nontrivial 
subspaces of Z and Z' are more exotic. For example, the Zaremba problem, where Dirichlet resp. Neumann 

conditions are imposed on two closed subsets Sjj resp. S^r of E with common boundary and covering E, leads 

_ i 

to the subspace V = K®H 2 (Ejv), which presents additional technical difficulties not covered by the ipdbo 
calculus. 

But when we go beyond the scalar second-order case, subspace situations have a primary interest; see the next 
section. 



3.3 Higher order operators and systems 

Let us now consider systems (matrix-formed operators) and higher order elliptic operators. Here one finds that 
subspace cases occur very naturally and allow studies within the V'dbo calculus, with much the same flavour as in 
Theorems 3.3, 3.4 and 3.6. For even-order operators, a general and useful framework was worked out in [17] — 
normal boundary conditions for operators acting between vector bundles — which could be the point of departure 
for ample generalisations. Scalar operators are covered by a simpler analysis in [16]. We shall here show in detail 
how the analysis of [16] can be used, and illustrate systems cases by examples. 

Example 3.11 Let A — (Aj^) be an p x p-matrix of second-order differential operators on fi, elliptic in the 
sense that the determinant of the principal symbol is nonzero for x 6 fi, £ ^ 0. We here have a Green's formula 
like (3.7), but where Sq(x) is a regular p x p-matrix and sq(x) is replaced by so(x)*. There is again a Dirichlet 
realisation Ay, and if it is elliptic and G q(Aj), we can repeat the study of Section 3.2, now for p-vectors. That 
will cover boundary conditions of the form (3.28), and give a somewhat abstract treatment of the more general 
cases. 

Now one can also consider boundary conditions where a Dirichlet condition is imposed on some components 
of Jqu and a Neumann condition is imposed on other components of v\u. This is very simple to explain when 
sq(x) = I, so let us consider that case. Let A be determined by a boundary condition of the type 

Iouno = 0, jiU Nl - C~f u Nl = 0, (3.71) 

where Nq and A^i are complementing subsets of the index set N = {1, 2, . . . ,p} with po resp. p\ elements. Using 
the homeomorphism 72 : Z ^ Ili<j<p (^0 an d similar vector valued versions of the other mappings in 
Section 3.2, one finds that realisations of boundary conditions (3.71) correspond to operators 

L: [] J] £T'(E), L = (7-pr JVi P 7 0i7li i JVll (3.72) 
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where pr^ projects {<fij}j£N onto {¥>j}jeJVii an d ijvi injects vectors {^>j)j^Ni into vectors indexed by N by 
supplying with zeroes at the places indexed by No. Again, ellipticity makes the realisation regular, and there are 
formulas very similar to those in Section 3.2: 

L X = L + WNA p L-yi ~ ^0,71)^ = C - P r "i P 7o,-n ^1. M ^ A ) = ~( lX )~\ ( 3 - 73 ) 
(A- A)- 1 = (A, - A)- 1 - K x i Nl M L (X)pi Nl (K'~ x )*, 

where L x and Ml (A) are elliptic ?/>do's when the boundary condition is elliptic. 

Notation gets a little more complicated if s is not in diagonal form, or if pr Wi is replaced by a projection onto 
a pi -dimensional subspace of C p that varies with x G S, Then it is useful to apply vector bundle notation, as 
in [17], where fully general boundary conditions are treated. 

Next, consider the case where A is a 2m-order operator A = J2\ a \<2m a a (x)D a (scalar or matrix-formed), 
elliptic on f2. The Cauchy data are the boundary values 

gu = {7 u, . . . , 7 2m -iu}, (3.74) 

which one can split into the Dirichlet data and the Neumann data 

JU = {70K, . . . ,7 m _iu}, uu = {^ m u,...,j2m-iu}. (3.75) 

There is a Green's formula for u, v G iJ 2m (f2): 

(Au,v) — (u,A'v) = (\u,jv) — (ju,x'v), x u = A\vu, x' v = A! x vv + AIq'JV, (3.76) 

where A\ and A[ (both invertible) and A' are suitable m x m matrices of differential operators over S. Here 
X = {xo, ■ ■ ■ , Xm-i] with Xj °f order 2m — j — 1, and the trace maps are continuous, 

j:Z s x (A)^ H ff s ^-5(S), II # s - 2m+i+ *(£), (3.77) 

0<j<m-l 0<j<m-l 

for s G M. If the Dirichlet problem is uniquely solvable, one can again use A 1 as reference operator, and set 

P x x = X K X ; P* x , = x'K; (3.78) 

with K x resp. K' x denoting the Poisson operator solving 

(A - X)u = resp. (A' -X)u = in fl, ■yu = if. (3.79) 

Now 

r = X - P^ x l and T'= X '- P'^l (3.80) 

are defined as continuous maps from D{A max ), resp. D(A' max ) to rij<m H^ + ^ (E), and there is a generalisation 
of (3.76) valid for all u G D(A max ), v G D(A' max ): 

(Au,v) - (u,A'v) = (ru,7«) { . + i }|{ _._i } - {ju,T'v){-j-i},{j+i} (3.81) 

(where (•, -){ Sj },{- Sj } indicates the duality pairing between Y[ J<m H s ' ( s ) and I\ J<m H ~ Sj 
It should be noted that P x x is a mixed-order system: 

P x x = (Pjk)j,k=o,...,m-r, Pjk of order 2m - 1 - j - k, 

and the principal symbol and possible ellipticity is defined accordingly. Such systems (with differential operator 
entries) were first considered by Douglis and Nirenberg [11] and by Volevich [44]; in the elliptic case they are 
called Douglis-Nirenberg elliptic. 
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Formula (3.81) can be turned into a boundary triplet formula with fi = K, = L 2 {Yj) m by composition of T 
and r' with a symmetric ?/>do defining the norms: 

6 = (^A fe+ i) fe , ;=0 m -i ■ 11 HJ+ H^) ^£ 2 (£) m , (3.82) 

0<j<m 

and composition of 7 with _1 , setting 

ri = er, r; = er', r = e _1 7 = r^ ; (3.83) 

this leads to a word-for-word generalisation of Proposition 3.1, with 

(Au,v) L2 (q) ~ {u,A'v) L2 (q) = (Tiu,r' v) L2 ( S)m - (ro^r'i^L^sjm, (3.84) 

foTueD(A max ),veD(A' max ). 

The straightforward continuation of what we did in Section 3.2 is the study of boundary conditions of the 
type x u = Cju, that are the ones obtained when A corresponds to T : Z — > Z' . When A corresponds to 
T : V — > W, we find that it represents a boundary condition similar to (3.58). Since the proofs in Section 3.2 
generalise immediately to this situation, we can state: 

Corollary 3.12 Theorems 3.3, 3.4 and 3.6 extend to the present situation for elliptic operators of order 2m, 
when H~ 2 (£) is replaced by Y[j< m ^ 3 ~ 5 (£), and 70 and v\ are replaced by 7 and X- 

Theorems 3.9 and 3.10 likewise extends when, moreover, ^(S) w replaced by L2(T,) m and (3.84) is used. 

Whereas the statements in the case of general subspaces V and W will be somewhat abstract, there are now 
also formulations where the subspaces are represented by products of Sobolev spaces, and the operators belong 
to the ?Adbo calculus. We shall demonstrate this on the basis of the treatment of normal boundary value problems 
in [16], which we now recall. Denote 

M = {0,...,2m-1}, M = {0,...,m- 1}, M ± = {m, . . . ,2m - 1}. 

A general normal boundary condition is given as 

f jU + B jklk u = 0, j G J, (3.85) 

where J is a subset of M with m elements and the Bj k are differential operators on E of order j — k. (It is called 
normal since the highest-order trace operators have coefficient 1 .) Let K = M \ J, and set 

J = JnM , Ji = JnMi, K = KnM , k 1 = kc\M 1 , 

l.h = {lj}j£Jo: V Ji = {7j}je./i, JK = {7j}j£K , v Kl = {-,},. K . 

then (3.85) can be reduced to the form 

7,7 it = F a ~f Ko u, v^u = Fij Ko u + F 2 v Kl u 1 (3.86) 

with suitable matrices of differential operators Fq, F\, F 2 . To reformulate this in terms of {7, \}, we use the 
convention for reflected sets: 

N' = {j I 2m - 1 - j e N}, (3.87) 

considered as an ordered subset of M. Then x = {Xj}jeM a splits into x = > XKi'}> where 

X.h' = {X]}]£Ji'i XKi' = {Xo)jeK 1 '- 
We can now reformulate (3.86) as 

j Jo u = Fq^k u, xji' u = GilK u + G 2 Xk 1 'U, (3.88) 
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where G\ and G2 are matrices of differential operators derived from those in (3.86) and the coefficients in Green's 
formula (more details in [16]). 

When A is the realisation of A determined by the boundary condition (3.88) (i.e., D(A) consists of the u £ 
D(A max ) satisfying (3.88)), and ellipticity holds, then D(A) C H 2m (Q), and the adjoint realisation A* is 
determined by the likewise elliptic boundary condition 

7kV v = -GZlJi'V, Xk v = Gil.h'v ~ F oX'j v- (3.89) 

Setting 

we have that 

1 D{A) cl = $( Yl H- k -?(Z)), jD(A*) cY = $( J] 
kei< jeJi' 

Here X is the graph of F and naturally homeomorphic to its "first component" Ylke k ^ k ~^ (^)> and similarly 
Y is homeomorphic to rijeji' H~ : >~i (£). The operator T : V — > that ^4 corresponds to by Theorem 2.1 
carries over to an operator L : X — > y * by use of 7, and this is further reduced to an operator 

Li ■ n ff_fe_ *(s) -> n = Gi - **p°, x $, 

fceffo je./i' 

with domain D(L\) = Y\ keKo H 2m ~ k ~^ (£). This representation is used in [16] to find criteria for the operator 
to be m-accretive, and the ideas are further pursued in [ 17] for systems of operators (where vector bundle notation 
is needed). 

For the present study of resolvents, we now find that A — A can be represented by 

L\ = Gi-**P 7 A X $ : " [ ff- fc -2(S)-> ] [ D{L\) = ] [ H 2m - k -^ (£), (3.90) 

when A G q(A 7 ), The corresponding A/-function and Krein formula are: 

(2 - A)" 1 = (Ay - A)" 1 - K^M(\)(K'^)*. (3.91) 

When A 6 f?(A) H g(A^), M(X) extends holomorphically to g(A) (note that the spectrum of £>(A 7 ) is discrete 
in this case). As a mixed-order operator with entries of order 2m — 1 — j — k (k £ A'o, j £ J\), L\ is 
Douglis-Nirenberg elliptic, and M(A) is so in the opposite direction. 

The two functions M(A) and A A together give a tool to analyse the spectral properties of A in terms of i/>do's 
on E, M(A) being holomorphic on and L A containing information on null-spaces and ranges. We have 
hereby obtained: 

Theorem 3.13 Let Abe a Im-order elliptic differential operator with coefficients in C°° (CI), and let A be the 
realisation defined by the normal boundary condition (3.85), reformulated as (3.88); assume that the boundary 
problem is elliptic. Let A 1 be the Dirichlet realisation, assumed elliptic and invertible. For A £ g(A 7 ), the 
operator corresponding to A — A by Corollary 2.3 carries over to (3.90). The associated M-function and Krein 
resolvent formula are described in (3.91); M(A) extends to an operator family holomorphic in A £ g(A). 

When A is scalar of order 2m, this analysis covers all the elliptic problems to which Seeley's resolvent 
construction [40] applies, for it is known that normality of the boundary condition is necessary for the parameter- 
ellipticity required there (cf. e.g. [21, Section 1.5]). When p x p-systems are considered, each line in (3.85) 
can moreover be composed with a multiplication map to the sections of a subbundle of £ x C p , as indicated in 
Example 3.11; here normality means surjectiveness of the coefficient of the highest normal derivative 7^, for each 
j. For such cases, the treatment can be based on [17]. 
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Example 3.14 For a simple example, consider the biharmonic operator A = A 2 , which has the Green's 
formula 

(A 2 «, v) - (u, A 2 v) = (xu, jv) - (ju, xv), xu = ( J 1 ^) ; (3.92) 



7 Am 

note that linear conditions on -foil, 71 u, 7 2 « and 73W can be written as conditions on 7oit, 71 u, 7oAu and 71 Au, 
hence on 71* and x u - The operator P A X is a t/jdo from H s ~^ x H s ~^ to H s ~^ x if s ~3 (note the reverse order 
of the Sobolev exponents in the target space). 

Let A be defined by an elliptic boundary condition with a first-order differential operator C, 

7 u = 0, 7 Au = C71M; (3.93) 

then A* represents the likewise elliptic boundary condition 

l0 u = 0, 7o Au = C*j lU , (3.94) 

and the domains of A and A* are contained in H 4 (fl). For the corresponding operator T : V — > W according to 
Theorem 2.1, V = W = K®({0} x H~*). This carries over to the operator 

L = C-pr 2 P° iX i 2 :H~i (3.95) 

where pr 2 = (0 l) and i 2 = ft^J . Thanks to the ellipticity, D(L) = Hi , and L is an elliptic ^>do of order 1. 
We here find that 

L x = L + pr 2 (P° - P*) i 2 = C- P r 2 i 3l Af L (A) = -(L A )-\ (3.96) 



7.X 7>X/ z ^ * 7>X " 

r 7 A i 2 M L (A)pr 2 (J, 7 



- A)- 1 = (Ay - A)- 1 - i 2 Mi (A) pr 2 (^ A ) 4 



It is also possible to take another realisation Ap than the Dirichlet realisation as reference operator, preferably 
one defined by an elliptic boundary condition, as in [15]. For strongly elliptic operators, using A 1 as reference 
operator has the advantage that semiboundedness properties are preserved in correspondences between A and T, 
see [16,17]. 

There do exist even-order operators for which the Dirichlet problem is not elliptic; a well-known example is 
the operator — A/ + 2 graddiv on subsets of IR 2 , studied by Bitsadze. 

Example 3.15 The operator A need not be of even order. For example, first-orderp x p-systems, such as Dirac 
operators, have received much attention. In this case the Cauchy data are {70U1, . . . , joUp}. Elliptic boundary 
conditions require p to be even. In "lucky" cases, one can get an elliptic boundary value problem by imposing the 
vanishing of half of the boundary values; this will then give a reference problem, allowing the discussion of other 
realisations. More systematically, one can impose the condition II + 7ou = for a certain ipdo projection II + 
over the boundary (the Atiyah-Patodi-Singer condition), which defines a Fredholm realisation. If it is invertible, 
it can be used as a reference operator. 

Similar considerations can be worked out for elliptic operators on suitable unbounded domains, and on man- 
ifolds. More specifically, the calculus of V'dbo's is extended in [21] to the manifolds called "admissible" there; 
they have finitely many unbounded ends with good control over coordinate diffeomorphisms. They include com- 
plements in R" as well as M. + of smooth bounded domains. 

Example 3.16 We shall here make concrete the considerations in Example 3.14 for the biharmonic operator 
on the half-space 

M-l = {x = {x',x n } 6 R n I x n > 0}, 

where we denote {xi, . . . , x n ~i} = x' . This constant-coefficient case can be viewed as a model for variable- 
coefficient cases, giving an example of the pointwise symbol calculations entering in the V'dbo theory. (A detailed 
introduction to the ijjdbo calculus is found e.g. in [22].) 
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It is well-known that since A 2 is symmetric nonnegative, the resolvent of the Dirichlet realisation in L 2 (R™ ) 
exists for A 6 C \ R+. We shall write 

A = n E K/4, where V e = {z E C \ {0} | | argz| < 9}. 

We fix a Ao < 0, then the realisations of A = A 2 — Ao fit into the general set-up, with the Dirichlet realisation 
A 1 of A 2 — Ao as reference operator. However, for simplicity in formulas (avoiding addition and subtraction of 
Ao), we keep the parameter A for the operator families defined relative to A 2 — A. 

In the following calculations, one can think of a /i > 0; the considerations extend holomorphically to V^/ 4 . 
To find the Poisson operator solving the problem 

(A 2 - X)u(x',x n ) = OonR'j, u(x', 0) = <p {x'), d n u(x', 0) = <pi(x'), 
we perform a Fourier transformation in the x'-variable, and then have to find solutions of 

((lef - dl) 2 + f i 4 )u(£,x n ) =0forx n > 0, 

d n 4(g,0) = M?), 0-91) 

that are in <S(R+) = r + iS(R) (the restriction to R+ of the Schwartz space of rapidly decreasing functions on R). 
Write 

am 2 - dif + m 4 ) = an 2 + *y - dimr - *m 2 - qd 

= (a + + d n )(a + - <9„)(<t_ + d n )(a_ - d n ), 

a + = (\a 2 + ^ 2 )K <J- = (\H'\ 2 -^ 2 )K (3.98) 

where is defined for z E to be real positive for z E R + ; note that /i 2 E V„/2 so that |£'| 2 ± i/i 2 E V^, and 
Re a± > 0. Then the general solution in <S(R+) of the first line in (3.97) is 

v(x n ) = Cl e~ a + Xn +c 2 e-°- x ". 

It is adapted to the boundary conditions by solution of 

ci + c 2 = <£o(£')> -a+ci - o-c-2 = 

with respect to (ci, C2); this gives 

Here k*(£' , x n ) is the so-called symbol-kernel of the Poisson operator Kij; it acts like 



andmapsff^^R™" 1 ) xiJ s -i(M™- 1 ) to {u E H S (R1) \ (A- X)u = 0} for alls E R. (Invariable-coefficient 
cases, k would moreover depend on x' .) 

I"n-Nf*iimnnn nnprafor } 

1,X 



The symbol p x (£,') of the Dirichlet-to-Neumann operator P x with \ chosen as in (3.92) is found by the 



calculation 



i^ 2 ( 2<T + <j- er+ + a- 



(7+ — <T- 



-(<T + +(T_) -2 

i/„ . \( <r+ + °-\_ (Poo Poi\ 



(3.99) 
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here we used that (cr + — er_) 1 = (cr + + cr_)/(o+ — cr?.) = (<r + + er_)/(4i/Li 2 ). 
The operator P A X is Op(p A ), where we use the notation for ipdo's 

Op(g)/ - (27T) 1 "" / e "' «'g(x', 0/(0 

Let us consider a model boundary condition (3.93) for A 2 — Ao with a first-order differential operator C = 
b±di + ■ ■ ■ + & n _i<9„_i with constant coefficients; it defines the realisation A, and C has symbol c(0 = ib ■ 
Then the operators L x ° and L x defined similarly to (3.95) and (3.96) are the ^do's with symbol 

/ A0 (O=c(O-Pn°(O, resp. 

l\0 = * A °(0 + P??(0 "Pn(0 = c(0 - Pn(0 = c(0 + + *_). 
When Z A is invertible for all Ml (A) is the -0do with symbol 

m(£', A) = -(c(0 + = -( C (0 + | a + V 2 )' + |(I0 2 -i** 2 )*)" 1 . (3-100) 

This shows how Ml (A) is found. As an additional observation, we remark that when b is real nonzero, then 
C* = -C ^ 0, so A is non-selfadjoint, cf. (3.94). However, Ml (A) is well-defined for all A g C \ 1+ 
(since c(0 G iR). On the other hand, L A ° has numerical range j/(L a °) approximately equal to the sector 
V = {z G C | | Imz| < |6| Re z}, and v(L x °) is contained in the numerical range of A, by [15, Th. Ill 4.3]. This 
gives an example where A has a large numerical range outside the spectrum. 

Remark 3.17 One can ask whether the considerations extend to Douglis-Nirenberg elliptic systems (systems 
of mixed order). But it may not be easy. Consider for example a 2 x 2-system 

on a smooth bounded open set f2 C R™, where u and v are scalar, and A+j is of order 4 — i — j. There is a Green's 
formula for A: 

(Aw,w') — (w,A*w') = (kw, 7ou') — {jou, k'w'), 

where kw = \u + bi^ov, k'w' = \'u' + &27ov', with first-order trace operators \ and x'\ {lo u i K w} are the 
reduced Cauchy data according to [13, 18]. Assume that the principal symbol is uniformly positive definite; in 
particular, the function A 2 2{x) > c > 0, the symbol a^x, £) > c|£| 2 and the determinant a\ x {x,£)A22{x) — 
ai 2 (x, a 2ii x > — C IC| 2 - Then the Dirichlet problem for A (with boundary condition jqu = 0) is well-posed, 
with domain D(A 1 ) = (H 2 {Q) n H%(Sl)) x H 1 ^). 

When G g(A~ ( ), there is the following parametrization of null-spaces of A (cf. also [18, Th. 2.1 1]): 

joWi ■ Z^iA) = {w G H S (Q) x H'-^Q^Aw = 0} ^ H s -i(Z), 
all s G R; but this does not in general lead to a nice parametrization of 

Z$(A) = {w G L 2 {VLf\Aw = 0}, 

which would be required to get a good interpretation of the abstract theory for A as an operator in L2(^) 2 - 
This regularity problem is not present in the one-dimensional situation, where the maximal domain is H 2 x H 1 ; 

an example is considered in Section 4.2. 

In the example in Section 4.1 we consider an n-dimensional case of (3.101), where the off-diagonal terms are 

of order 0; this allows an easier parametrization of the null-space. 



Copyright line will be provided by the publisher 



mn header will be provided by the publisher 



27 



4 Examples 

We have seen that the family Mj(A) is holomorphic on g(A) so that 

a{A) D {A e C | M A {\) singular at A}. (4.1) 

Among the singular points, we have very good control of those outside g{Ap) by Corollary 2.12 and (2.43), for 
the null-space and range of A — A are fully clarified by the same concepts for T A (a holomorphic operator family 
on g(Ap) which is homeomorphic to the inverse of Mj(A) on g(A) n g{Ap)); it gives the information: 

a{A) n g(A ) = {A e g{A f3 ) | kcrT A ^ {0} or ranT A ^ W- x }. (4.2) 

So the only spectral points of A whose spectral nature may not be controlled by Mjj(A) and T A are those that 
lie in a(Ap). For many scalar equations it has long been known that the M-function allows full control of the 
spectrum. However, when considering systems, it is easy to see that uncontrolled points may exist by considering 
equations which are decoupled. In [7], an example involving a coupled system of ODEs was given where 
was regular at a point Ao belonging to the essential spectrum of A (and of Ap). We shall here show a similar 
phenomenon for PDEs and for a system of ODEs with first order off-diagonal entries. 



4.1 PDE counterexample 

Consider the 2x2 matrix-formed operator 

A a(x) 



b(x) c(x) J ' (43) 

acting on 2-vector functions w = {it, v} on Q, such that Aq is a second-order elliptic operator as studied in 
Section 3.2 and a,b,c £ C°°(f2). The set ranc is a compact subset of C. We assume that it has a connected 
component K with more than one point, that C \ (ran c) is connected, and that 

a(x)b(x) vanishes on supp c. (4.4) 

The maximal and minimal operators are 

..in ax a(x)\ A ( A m in a(x) \ , . 

b(x) c{x) ) ' min ~ ^ b(x) c(x) J ' ( ' 

with 

D(A max ) = D{A , max ) x L a (fi), D(A min ) = D{A , min ) x L 2 (0), (4.6) 
and there is the Green's formula 



u\ I u \ \ I I u 
v ' 



with notation as in (3.7) ff. 

Proposition 4.1 Let A^ be the Dirichlet realisation defined by the Dirichlet condition 70 w = 0. 

(i) A~f is lower semibounded with domain (H 2 (fl) R Hq(CI)) x L2(fi), and the spectrum is contained in a 
half-space {Re A > a}. 

(ii) The operator A 7 has a non-empty essential spectrum, namely 

(7cs S (A 7 ) = ranc. (4.7) 

(iii) Outside the essential spectrum, the spectrum is discrete, consisting of eigenvalues of finite multiplicity. 
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Proof, (i) follows from standard results for elliptic operators and the fact that the adjoint of A 1 is the Dirich- 
let realisation of A! with similar properties. Hence, the spectrum is contained in a half-space {Rc A > a}. 

(ii) follows from Geymonat-Grubb [13], where it was shown that the essential spectrum of the realisation of a 
mixed-order system A defined by a differential boundary condition f3u = consists exactly of the points A £ C 
where ellipticity of {A — A, (3} fails. The current operator A is a mixed-order system with orders 2 and for the 
diagonal terms, order 1 for the off-diagonal terms (to fit with the rules of Douglis, Nirenberg and Volevich for 
mixed-order systems). The principal symbol of A — A is 

a°(*,0 (4g) 



c(x) - A 

so ellipticity of the Dirichlet problem for A — A fails precisely when A £ ran c. 

For (iii), we note that the resolvent set is non-empty, so for A ^ ran c, A 1 — A is a Fredholm operator (by the 
ellipticity) with index (since the index depends continuously on A). Then all spectral points outside ranc are 
eigenvalues with finite dimensional eigenspaces. Since these eigenspaces are linearly independent, there can only 
be countably many, so there is at most a countable set of eigenvalues outside ran c. They can only accumulate at 
points of ran c. □ 

Consider another boundary condition for A, 

v x u = C70U, (4.9) 

with C a first-order differential operator, and such that the system {Aq, v\ — C70} is elliptic, defining the real- 
isation Aq. Then {A, (v\ — Cjq) W\\ is likewise elliptic, and we define A to be the realisation of A under the 
boundary condition (4.9). Again, the essential spectrum equals ranc. 

It is well-known that Aq satisfies a 1-coerciveness inequality (hence is lower bounded) if and only if the real 
part of the principal symbol of L = C — P® Vl is > co|£'| with cq > 0, and that the adjoint then has similar 
properties (cf. e.g. [16], [17]). Assuming this, we have that the spectrum lies in a half-plane {Re A > ai}, and 
the spectrum is discrete outside ran c. 

We next want to discuss A/-functions for the comparison of the Dirichlet realisation and the realisation defined 
by (4.9). Let k be a point in g(A 7 ); then the general analysis of Section 2 works for the realisations of A — k, 
with reference operator A 1 — k. So the holomorphic families T x and M,£_ k s (A) are well-defined relative to 

this set-up. Note that c(x) is now replaced by Ck(x) = c(x) — k; the essential spectra of A 1 — k and A — k are 
contained in rancj, = ranc — k. 

Let A ^ rancfc. The solutions of the Dirichlet problem with non-homogeneous boundary condition are the 
solutions of 

(Aq — k — X)u + av = 0, , 
bu + (c — k — X)v = 0, 

7o" = <P- (4.10) 

The second line is solved by v = b(X + k — c) _1 w, which by insertion in the first line gives 

(A - k - A + ab(\ + k- c) _1 )?i = 

j Q u = ip. (4.11) 

When A £ g(A 1 — k), the problem (4. 11) has a unique solution u — Ky X tp £ ker(ylo~ k — \+ab(\+k— c) _1 ) 
for each 93 £ H~^, and (4.10) has the solution {u, b(k + A — c) _1 w} £ ker(A 1 — k — A). The Dirichlet-to- 
Neumann operator family for A — k is 

P k > x = v lWl K*> x , (4.12) 

which identifies with the Dirichlet-to-Neumann operator family Ph Ul for Aq — k. 
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As above, let A be the realisation of A under the boundary condition (4.9). It corresponds as in Section 3.2 to 

L x = C - P KX = L + P kfi - P k - X : H~i -> H%, where 
L = C-P k -°, D{L) = Hi. 
Then there is an M-function going in the opposite direction and satisfying 

Mi (A) = —{L + P fc '° - P^)" 1 (4.13) 
when A £ g(A — k) n f?(Ay — fc). With a notation similar to (3.14), Mi (A) acts as follows: 

Mi(A) = P^ C70i7Q . (4.14) 

Proposition 4.2 There exists a point belonging to the essential spectrum of A — k (and of A^ — k) at which 
Mi (A) is holomorphic. 

Proof. Let Ao = c(xrj) — k, for some xo where cq(x) £ K \ {0}. Then Ao belongs to a css (A 1 — k) and to 
a ess (A — k). We shall show that Mi (A) can be extended holomorphically across Ao or a neighbouring point. 
In view of (4.4), 

a{x)b{x){\ + k - cix))- 1 = a(x)b(x)(X + k)~ 1 , for all x e 0, all A ^ -k. (4.15) 
This implies that the problem (4. 1 1) takes the form 

(A - k - A + ab(\ + fc) _1 )u = 

7o" = <P, (4.16) 
for A 7^ — k, and this is obtained by reduction from the problem 

(Ao — k — X)u + aw = 0, 
bu - (A + fc)u = 0, 

70" = ^, (4.17) 

a Dirichlet problem for 

A 1 =( A °~ k a(X A- (4-18) 



b(x) —k 
Similarly, the problem 

(A — k — X)w = 0, i/iu - C70W = V, (4.19) 
is equivalent with 

(A t - X)w = 0, vxu - C70U = ip. (4.20) 
So, Mi (A) defined above coincides with the analogous operator for A%: 

Mi(A) = P^- Cj0 ,-y (A^ ( 4 - 21 ) 

It is holomorphic on q(A\ ), where A\ is the realisation of A\ defined by the boundary condition (4.9). 

This detour via A\ gives information about the possible holomorphic extensions of the Mi(A)-function for 
A — k. We infer from the general result of [13] that the Dirichlet realisation A\ tl of A\, as well as the realisation 
A\, have essential spectra equal to { — k}. Moreover, their spectra are contained in a half-space, and are discrete 
outside the point {— k}. 

So Ao = c(xq) — k is either in q(A\) or is one of the discrete eigenvalues of A\, and in any case there is a disk 
B(Xq, 5) around it such that Ml(A) is holomorphic on B(Xo, S) \ {Ao}. Since c(xo) — k is not the only point in 
the connected set K — k, there will be a point x' Q such that c(x' ) — k £ B(Xq, 5) \ {Ao}. 

We can conclude that Mi (A) is holomorphic at X' = c(x' ) — k, but the point belongs to the essential spectrum 
of A - k (and of A 1 - k). □ 
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Remark 4.3 The hypothesis (4.4) can be replaced by a weaker hypothesis as follows: Assume that K is 
a compact connected subset of ranc containing more than one point. Let uj,uj',uj" be subsets of il with 
dist(w,n \ uj') > 0, dist(a/,n \ uj") > 0, such that K C c(ui) and dist(K, c(Ti \ uj')) > 0. Assume that 
ab vanishes on uj" . Let rj G C°° with r\ = 1 on ft \ w" and 77 = on a/, and set c' = 77c. Then, for all 
A ranc' — k, ab(\ + k — c) _1 = ab(X + k — c') _1 . The problems (4.10) and (4.19) can now be replaced by 
problems where c is replaced by c' , whose range is disjoint from K, so that there will be points in K — k where 
the ill-function is holomorphic. 



4.2 ODE counterexample 

Consider the 2x2 matrix-formed operator 



(x) 



—u"{x) a(x)v'(x) 
b(x)u'(x) c(x)v(x) 



(4.22) 



acting on pairs of functions u, v on the interval [0, 1] and a, b, c G C°°([0, 1]). Its formal adjoint is given by 

-u"(x) (b(xjv(x))' 



A' 



Or) 



^(a(x)u(x))' c(x)v{x) 
and there is the Green's formula 



(4.23) 



where 



.4 



Fiit 



(r lU , r' v) - (r «, i» 



-ui(l) + a(l)u2(l) 
ui(0) -a(0)u 2 (0) 



, T u = T' u 



ui(l) 
ui(0) 



r> 



-«i(l)-6(l)«a(l) 
«i(0) + 6(0)^(0) 



(4.24) 



andii,7j g D(A max ) = H 2 {0,1) x H 1 ^,!). 
Proposition 4.4 We introduce the quantities 



Q(x,X) 



-1 



a(x)b(x) 
\ — c(x) 



, a(x) = exp / Q(s, A) ds 1 , /3(a;) 



a(x) 



-1 



a(a:)6(j:) 
A — c(:r) 



anc/ /ef y\ and 1/2 be two linearly independent solutions of 
1 

4a 2 



7? 



/3A 



= 



(4.25) 



(4.26) 



satisfying the initial conditions ?/i(0) = l,j/^(0) = 0, 2/2(0) = 0,y 2 (0) = !■ Note that yi and y 2 depend on A, 
£>wf we suppress this in the notation. 
(i) We /?ave 



ker(A n 



A) 



cij/i + c 2 y 2 
rMciyJ + c 2 y 2 ) 



ci,c 2 G 



(ii) Consider the operator A\, the restriction o/^4 max to ker Y\, i.e. subject to Neumann boundary conditions. 
The matrix of the M -function has the entries 

2/i(l) 1 



mn 



m 2 i 



a(i)fa(i) 

A-c(l) 



- 1 



2/1(1) 



-, mi2 



1 



"(i)fe(i) 

A-c(l) 



1 



TO22 



a(0)o(0) 
A-c(0) 

2/2(1) 



2/1(1) 



1 



a(0)b(0) 
A-c(O) 



2/1(1) 



(4.27) 



(iii) Assume in addition that the function b vanishes identically on an open interval I C (0, 1). Then there is a 
second operator A with different essential spectrum to A\, but giving rise to the same M -function. 
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Remark 4.5 Note that if y[ (1) = 0, then yi is a Neumann eigenfunction for (4.26), and A is an eigenvalue of 
the operator A%. Apart from this case, singularities of the M -function only occur when the coefficient in (4.26) 
blows up, as terms of the form a ^_^^ — 1 also appear there. 

Proof, (i) (A — X)u = can be written as 

— u'{ — Xu\ + au 2 = 0, bu 1 + (c — A)« 2 = 0. 
Solving the second equation for u 2 and substituting this into the first gives 

bu\ y 

—j- J - Xui = 0. (4.28) 
Introducing a and 3 as in (4.25), the equation (4.28) simply becomes (au^)' — BXui = 0. Moreover, introducing 

/x 
— U!(x), 

the equation can be written as 
1 



a 



a' 2 a" ox 

BX 

4a 2 1 



<p = (4.29) 



and the kernel of A max — A has the form 

C12/1 + C 2 ?/2 



ker(A max - A) 



I^(ci2/i +c 2 y' 2 ) 



ci,c 2 G 



(ii) We now calculate the M-function for the operator subject to Neumann boundary conditions. For any 

(ui, u 2 ) S ker(A max - A) we have 

fm n m 12 \ f-u[{l) + a(l)u 2 (l)\ = (ui{l) 
\m 2l m 22 ) \ ui(0)-a(0)u 2 (0) j \u x {Q) 

A simple calculation then gives (4.27). 

(iii) We first note that the determinant of the principal symbol of A — X is given by £ 2 (ab + c — A) which is zero 
if A S ran(a6 + c). Hence ellipticity of the system fails for A £ ran(a6 + c) and by [13], the essential spectrum 
of the operator equals ran(a6 + c). Under the additional assumption that the function b vanishes identically on 
an open interval / C (0, 1), we have that ran(c| /) is contained in the essential spectrum of the operator. Let c be 
a C°°-function on [0, 1] that coincides with c on [0, 1] \ /. This gives rise to another operator A with 

A (*) = (uSl tl V lS) (4-30) 



v J \b(x)u'(x) c(x)v(x) 

Let A\ be the realisation of A subject to Neumann boundary conditions. Then ran(cj/) will lie in <r css (Ai). 
Therefore, in general the essential spectrum of the operators Ai and A\ will differ. However, for the calculation 
of the M -function, c only appears in terms of the form ^^(x) an< ^ our assumptions, 

b(x) b(x) 



X — c(x) X — c(x) ' 

so the M-functions for A\ and A\ coincide. Thus, we have another example where two operators with differing 
essential spectra give rise to the same A/-function. □ 
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